Nonequilibrium Dynamics of Anisotropic Large Spins in the Kondo Regime: 
Time-Dependent Numerical Renormalization Group Analysis 



00 
O 

o 



David Rooscn, 1 Maarten R. Wegewijs, 2,3 and Walter Hofstetter 1 

1 Institut fur Theoretische Physik, Johann Wolfgang Goethe- Universitat, 60438 Frankfurt/Main, Germany 
z Institut fur Festkorper-Forschung - Theorie 3, Forschungszentrum Jiilich, 52425 Jiilich, Germany 
3 Institut fur Theoretische Physik A, RWTH Aachen, 52056 Aachen, Germany 

(Dated: May 25, 2007) 

We investigate the time- dependent Kondo effect in a single-molecule magnet (SMM) strongly 
coupled to metallic electrodes. Describing the SMM by a Kondo model with large spin S > 1/2, 
we analyze the underscreening of the local moment and the effect of anisotropy terms on the re- 
laxation dynamics of the magnetization. Underscreening by single-channel Kondo processes leads 
to a logarithmically slow relaxation, while finite uniaxial anisotropy causes a saturation of the 
SMM's magnetization. Additional transverse anisotropy terms induce quantum spin tunneling and 
a pseudospin-1/2 Kondo effect sensitive to the spin parity. 
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Introduction. Crystals of single-molecule magnets 
(SMMs) are mcsoscopic systems, which exhibit funda- 
mental quantum mechanical phenomena like quantum 
tunneling of the magnetization (QTM). The individual 
SMMs responsible for these effects have large spin S ~ 10 
and show rich quantum dynamics due to the anisotropic 
effective potential induced by ligand-fields and spin-orbit 
coupling [1]. More recently, transport through individ- 
ual SMMs coupled to metallic electrodes has been mea- 
sured for the first time [2, 3], giving access to many-body 
physics in these complex quantum impurity systems. In 
particular, the large spin, in combination with magnetic 
anisotropy, was predicted to give rise to a pseudospin- 
1/2 Kondo effect [4] in a single SMM. In combination 
with applied magnetic fields, this effect can be used for 
transport spectroscopy of an individual SMM in a two- 
terminal geometry [5] . Beyond the linear response trans- 
port investigated in these works, one can gain fundamen- 
tal insight by studying nonequilibrium magnetization dy- 
namics of an SMM coupled to macroscopic electrodes. 
Recent STM experiments have also demonstrated that 
single molecules attached to a substrate can be abruptly 
switched to a new conformation with different transport 
properties by a voltage pulse, resulting, for instance, in 
a change in the Kondo temperature [6] . Recent advances 
in the field of quantum computing using electron spins 
in semiconductor quantum dots concern high-frequency 
control of magnetic and electric signals. In particular, 
recently it has been shown experimentally [7] that via os- 
cillating electric fields spin transitions can be induced lo- 
cally by spin-orbit coupling. There are prospects for ap- 
plying such techniques to three-terminal single-molecule 
devices although technical challenges remain. 

Motivated by these advances, here we theoretically 
analyze the nonequilibrium spin dynamics of a SMM 
in response to a sudden change of the magnetic field 
strength, with a special focus on the manifestation of 
Kondo screening in the time domain. The accurate real- 




FIG. 1: Energy diagram of the SMM. The levels indicate spin 
eigenstates \S,M) which are split in energy by the uniaxial 
(D) anisotropy. Two types of spin fluctuations occur: pro- 
cesses due to the intrinsic transverse anisotropy (B2) , result- 
ing in QTM, and exchange processes involving the conduction 
band electrons (J), leading to the Kondo effect. 



time calculation of impurity observables like the magne- 
tization is crucial for developing a general understand- 
ing of time-resolved Kondo spectroscopy of SMMs, for 
instance using transport measurements [5]. More gen- 
erally, proposals for spintronic operation of SMM-bascd 
devices [8, 9, 10, 11, 12] rely on a detailed understand- 
ing of the magnetization dynamics. Note that here we 
consider pulsed magnetic fields as opposed to the exper- 
imentally more easily accessible oscillating fields, since 
the former are at the moment more conveniently treated 
using numerical renormalization group (NRG) methods 
[13, 14, 15, 16]. In order to analyze the effects of pulsed 
fields we use the time-dependent numerical renormaliza- 
tion group method (TD-NRG) recently introduced by 
Anders and Schiller [16]. The advantage of this approach 
is that the time evolution of any quantum impurity sys- 
tem after a sudden perturbation at t = can be calcu- 
lated in a numerically exact way on arbitrarily long time 
scales without accumulating any discretization error that 
scales with the elapsed time, in contrast, e.g. to the 
adaptive time-dependent DMRG [17]. For a detailed de- 
scription of the TD-NRG method we refer to [16]. For all 
the results presented below we have exhaustively checked 
that the dependence on the method's parameters is negli- 
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gible (such as the discretization parameter A or the num- 
ber of states NL eve i taken into account in each iteration) . 
We employ units in which g = H = fcg = W = 1, where 
W is half the width of the conduction band. 

Model. We consider a SMM as a (generalized) quan- 
tum impurity coupled to a conduction band of metallic 
electrodes : H = H c b + 'Hint + T~tim P [4]. The conduc- 
tion band and the interaction term in the Hamiltonian 
read H c b = Hk,u e k c \ u c kfl ,TL lnt = J S ■ s, where s is 
the spin of the most localized electronic orbital in the 
conduction band which is coupled to the impurity spin 
S (e.g. S = 10 for the single-molecule magnet Mni 2 
and S = 8 in case of Fes). Below we systematically in- 
vestigate different models describing the impurity 7ii m p 
(see Fig. 1). The time scale on which physical quantities 
vary throughout this work typically ranges from 10 3 to 
10 9 times h/W corresponding to 1 ps - 1 [is for typical 
values of the spin, exchange interaction and anisotropy 
parameters. 

Isotropic Kondo model. Let us first investigate an 
isotropic spin S with antiferromagnetic exchange J > 
induced by electron tunneling, i.e. Tti mp = 0. For a 
single electron with spin S = 1/2 at zero temperature 
and zero magnetic field the local moment is completely 
screened by the Kondo effect, i.e. (S z ) = 0. Here, how- 
ever, we consider a large impurity spin 5* > 1/2, as rel- 
evant for SMMs [4], aligned at times t < by a strong 
magnetic field h z in the longitudinal direction, which is 
sufficiently strong to completely polarize the impurity 
spin. At t = the field is switched off and using the 
TD-NRG we calculate the screening of the SMM's mag- 
netization by exchange scattering processes as a function 
of time, see Fig. 2. The most prominent observation is 
that, in contrast to the S = 1/2 case, the interaction 
between the impurity and the conduction electrons does 
not quench the magnetic moment on the impurity com- 
pletely. This is referred to as the underscreened Kondo 
effect [18, 19, 20, 21, 22]. In a naive picture the equi- 
librium underscreened Kondo effect can be described as 
a two-stage process: first the magnetization is reduced 
to S — 1/2 by screening due to the spin-1/2 electrons in 
the conduction band. This effective spin S' = S — 1/2 
then couples ferromagnetically (</' < 0) to the conduc- 
tion band and becomes asymptotically free at low en- 
ergies [19]. Our results in Fig. 2 show how this energy 
dependence of the effective couplings translates into the 
real-time dynamics of the impurity spin. For different 
values of S and J = 0.4 we observe an initial drop of 
the magnetization by approximately 1/2, followed by a 
logarithmically slow further relaxation (note that (S z ) is 
not a conserved quantity and can drop below the value 
of S z — 1/2). Comparing the low energy excitations for a 
spin S and antiferromagnetic coupling to those of a spin 
S' = S — 1/2 impurity coupled ferromagnetically, Kollcr 
et al. [22] have determined the value of the renormal- 
ized ferromagnetic interaction J' for the underscreened 
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FIG. 2: The logarithmic decay of (S z )(t) due to underscreen- 
ing of the local moment for different spin sizes S (J = 0.4) 
is compared to the relaxation of a spin S' = S — 1/2, whose 
effective ferromagnetic coupling J' < is chosen according to 
[22]. (NRG parameters: A = 2.5, N = 60, N Leve i = 2000). 

Kondo effect. In Fig. 2 we show that this correspon- 
dence can also be observed in the real-time dynamics, 
by comparing the underscreened J = 0.4 relaxation dy- 
namics with a ferromagnetic S' = S — 1/2 Kondo im- 
purity and an effective J' < as calculated in [22]. In 
the ferromagnetic case, we fitted the time-evolution as 
(S z (t)) = A/[log(t + B)] using A and B as fit parameters 
and found very good agreement with our results. For 
antiferromagnetic coupling good asymptotic agreement 
with the corresponding ferromagnetic case was found in 
the time evolution. This confirms our identification of 
time-resolved underscreening. We observe that the time 
scale of the partial screening 5^5 — 1/2 increases with 
the spin size S. If the exchange interaction J is decreased, 
the effect is more dramatic: the partial screening sets in 
at larger times and its time scale increases considerably 
(not shown) in agreement with the non-linear increase of 
the effective coupling J' for small J as found in [22]. 

Kondo model with uniaxial anisotropy. Single-mo- 
lecule magnets are characterized by additional magnetic 
anisotropy terms in the impurity part of the Hamilto- 
nian, as described by the generic model TLi m p = ~DS Z — 
\B 2 [S\ + S^_) — h ■ S. The spin has to overcome an en- 
ergy barrier of height DS Z to reverse its ground state 
alignment with the easy axis intrinsic to the SMM. 
Ligand field effects reduce the symmetry group of the 
molecule about this axis down to a finite group. The 
non-commuting terms oc B2 allowed by this low sym- 
metry introduce quantum tunneling of the magnetiza- 
tion (QTM) through the barrier [1]. The energy levels of 
an isolated SMM are sketched in Fig. 1. The resulting 
model has recently been studied [4, 5] since it is relevant 
to transport experiments on SMMs immobilized between 
metallic electrodes [2, 3]. 

We first investigate how the development of the under- 
screened Kondo effect is hampered by a finite anisotropy 
barrier D > in the absence of QTM (B 2 = 0). When 
switching off the magnetic field, a clear saturation of 
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FIG. 3: Suppression of underscreening: the magnetic 
anisotropy barrier D > impedes further relaxation of the 
SMM's magnetization by the Kondo effect (S = 3/2, J = 0.6). 
(NRG parameters: A = 1.8, N = 60, N Levet = 1000.) 
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the magnetization on a time scale oc 1/D is seen in 
Fig. 3. The logarithmic decay characteristic of the un- 
derscreened Kondo effect without anisotropy D is termi- 
nated, because Kondo spin fluctuations with small energy 
uncertainty are suppressed by the barrier. In a renormal- 
ization group language, the scaling is cut off at D and the 
Kondo screening stays incomplete. 

Kondo model with uniaxial and transverse anisotropy. 
The transverse anisotropy term B 2 induces QTM in 
SMMs in the absence of a magnetic field. We first 
consider the typical case B 2 <C D and strong coupling 
to the conduction band J = 0.2 for half-integer spins 
S > 1/2. The Kondo effect developing in this limit 
was shown [4] to involve all magnetic excitations of the 
SMM, since Tk{S, J, B 2 , D) exceeds the magnetic split- 
tings due to the large J value. As before, we prepare a 
spin 5 = 3/2 system in an almost completely polarized 
state {S z )(0) « 3/2 at t < 0, and monitor the subsequent 
time dependence of the magnetization (S z )(t). For very 
weak transverse anisotropy (D = 5 ■ 1CP 7 , B 2 = 8 ■ 10~ 9 , 
J = 0.2), three different regimes can be distinguished in 
the time evolution (cf. Fig. 4): (i) Starting at t rts 1/W 
the conduction band electrons partially screen the spin on 
the impurity. Both the time evolution and the NRG en- 
ergy levels are the same as for the "bare" underscreened 
Kondo effect obtained by setting B 2 = D = 0. (ii) 
The transverse anisotropy B 2 introduces a new feature: 
for times tW > 100 (i. e. before the finite barrier D 
"freezes" the spin) the magnetization starts declining and 
displays an oscillation around the value (S z ) ~ 0.35. (iii) 
The remaining magnetization is quenched by spin ex- 
change processes between the conduction band electrons 
and the impurity. The time scale is given by the inverse 
Kondo temperature of the pseudospin-1/2 Kondo effect 
r oc 1/T K (S, J, B 2 ,D) in line with [4, 5]. The small rem- 
nant magnetization at large times is due to incomplete 
thcrmalization which is unavoidable in a finite size sys- 
tem with discretization parameter A > 1; see also [16]. 



FIG. 4: Half-integer spin (J = 0.2). Three different regimes 
in the time evolution relate to underscreening (tW < 100), 
buildup of the pseudospin-1/2 Kondo effect (100 < tW < 
10 6 ), and approach to the strong-coupling fixed point (tW > 
10 6 ). Above the time evolution we plot the eigenenergies of 
the NRG Hamiltonian H m as a function of iteration m, i.e. 
the spectrum of the many-body system SMM + electrode, 
where conduction-band electrons with energies above 1/A m 
have been integrated out exactly. The scales for the two plots 
are adjusted to each other, allowing a direct comparison be- 
tween the NRG level flow and the corresponding real-time 
evolution. (NRG parameters: A = 2.0, A" = 70.) 



In order to study the sensitivity of the observed mag- 
netization dynamics to the spin size S as well as the ex- 
change coupling </, we now consider large values for the 
anisotropy barrier D ~ lO -2 !^ 7 and the quantum tunnel- 
ing term B 2 ~ D, for a weakly coupled SMM (J = 0.1). 
For half-integer spin S the quantum tunneling should 
then reduce the magnetization to a finite value, before 
the underscreened Kondo effect sets in. Our results for 
the time-dependent magnetization (S z }(t) indeed display 
this behavior: Fig. 5 (a) shows a clear two-step relaxation 
and a strong dependence on the half-integer value of S 
as well as on the exchange coupling J. At tW « 10 the 
quantum tunneling term B 2 starts mixing the different 
magnetic states of the molecule. The magnetization is 
reduced to (S z ) « 0.7 showing damped Rabi-typc oscil- 
lations around this value for all spin sizes we investigated 
(S = 3/2 — 7/2), which are related to the quantum tun- 
neling term and which disappear for B 2 = 0. For larger 
times t, we again find the pseudospin-1/2 Kondo effect 
to set the relevant time scale r oc 1/Tk{S, J, B 2 , D) for 
the complete screening of the molecule's magnetization. 
Since the Kondo temperature increases with the size of 
the spin (cf. [4]) the relaxation is faster for larger spin 
values S [cf. Fig. 5 (a)]. The dependence of the relax- 
ation on the exchange coupling J for fixed spin S = 5/2 
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FIG. 5: Magnetization dynamics for weak Kondo coupling 
J = 0.1, large QTM, B 2 = 0.75 • D = 7.5 ■ 10~ 3 , and different 
sizes of the spin: (a) Half-integer spin: two-stage relaxation, 
with a time scale for the Kondo effect strongly dependent on 
the size of the spin. Inset: dependence on exchange coupling 
J [shown for spin S = 5/2, J = 0.001 (red) and J = 0.05 (or- 
ange) ...0.15 (purple) in equidistant steps], (b) Integer spin: 
the different spin sizes have only little effect on the damped 
Rabi oscillation of the magnetization around zero. (NRG pa- 
rameters: A = 2.0, N = 65.) 



is shown in the inset of Fig. 5 (a). For larger J the oscil- 
lations reach negative values. The narrow plateau after 
the first "drop" in the magnetization, which is due to 
the transverse anisotropy £?2, saturates at (S z ) « 1.5, 
and the time scale for the pseudospin-1/2 Kondo effect 
becomes very large in the limit of weak coupling J. 

For integer spin values S the QTM term alone al- 
lows the magnetization to tunnel between the two ground 
states IS") and \-S) of the isolated SMM (J = 0), result- 
ing in zero net equilibrium magnetization at low magnetic 
field. Therefore, one does not expect a Kondo effect, be- 
cause no cooperation of the spin exchange processes is 
needed to change the direction of the magnetization [4] . 
We indeed observe a dramatic change in the spin dynam- 
ics going from half-integer spin to the nearest full-integer 
S for identical parameters B2, D and J [compare Fig. 5 
(a) and (b)]. {S z )(t) displays an oscillation around zero 
which dies out on time scales ~ 10 4 /!! 7 due to damping 
induced by the conduction band electrons. 



Conclusions. In this work we have analyzed the 
noncquilibrium magnetization dynamics of Kondo mod- 
els with large spin S > 1/2 in response to a sud- 
den perturbation. For the isotropic Kondo model we 
found a two-stage relaxation, which displays the under- 
screened Kondo effect in the real-time magnetization dy- 
namics and confirms that the renormalized coupling to 
the impurity after a partial screening of the magneti- 
zation to S — 1/2 is ferromagnetic J' < 0. For the 
extended anisotropic Kondo Hamiltonian describing a 
single-molecule magnet, wc found a suppression of the 
underscreened Kondo effect by uniaxial anisotropy. With 
additional transverse anisotropy, relaxation becomes sen- 
sitive to the spin parity. For half-integer spin, a complex 
interplay of tunneling and spin-screening by conduction 
electrons leads to a pseudospin-1/2 Kondo effect and a 
reduction of the magnetization on long time scales. In 
contrast, for integer spin, the transverse anisotropy only 
leads to damped Rabi oscillations without spin screening. 

The authors thank Christian Romcikc for discussions. 
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